Decoherence of a Single Qubit from Quantum Anomaly 
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We show via an explicit example that quantum mechanical anomalies can lead to decoherence 
, of a single quantum qubit through phase relaxation. The anomaly causes the Hamiltonian to 

■ develop an non-self-adjoint piece due to the non-invariance of the domain of the Hamiltonian 

04 ' under symmetry transformation. The resulting decoherence originates completely from the 

dynamics of the system itself and not from interactions with the environment. 

(N 



A physically realizable quantum computer must satisfy some dehcate requirements One of 
these requirements is that coherence must be maintained within a single qubit and also among 
entangled qubits. Coherence within a single qubit requires dynamics of the two-level quantum 
state to be controlled by unitary evolution. This in turn is guaranteed by the self-adjointness of 
the Hamiltonian in the Schrodinger equation. Up to now, attention has been given mostly to deco- 



^r~i . herence that originates from the interaction of the quantum system with its external environment 

' 2]- The purpose of this letter is to point out that decoherenece can also come from anomalous 

, symmetry breaking of the quantum mechanical system. The novelty of this phenomenon is that 



the decoherence originates from the system itself and not via interactions do with an external en- 
vironment. This anomalous decoherence, which we make explicit in the following via a toy model, 
is potentially significant for quantum information theory and should in principle be taken into 
account in the construction of quantum computing models. 
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H The model we consider is an electron in a magnetic field produced by the Dirac monopole. It 

is described by the following Hamiltonian, 

2m 2m 2m 

where A is the singularity-free vector potential of the Dirac magnetic monopole P] and B = 
V X A = gr/r^, the corresponding magnetic field. Further, the single- valueness of the wave 
function requires that 2eg should be an integer. 

The model ([TJ possess a rotational symmetry SO (3) as well as a dynamical superconformal 
symmetry OSP{l,l) [4]. The S0{3) symmetry is generated by the angular momentum of the 
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electron-monopole system, J = r x (p — eA) — egr/r + a /2. The 05*^(1, 1) consists of two parts: 
one is the conformal symmetry S0{2, 1) generated by the Hamiltonian H, the dilatation operator 
D and the conformal generator K, and the other part is the N = 1/2 conformal supersymmetry 
generated by the supercharge Q and conformal supersymmetry generator S. The S0{2, 1) con- 
formal symmetry is a generic feature of physical systems with potential [S], whose algebra is 
reahzed as [H, D] = iH, [H, K] = 2iD, [D, K] = iK. 

The large symmetry described above allows the model ([1]) to be solved exactly with a suitable 
representation of S'0(3) x 0S'P(1, 1) |4i]. The quantum states are characterized by the eigen- 
states \j,m,a,E) of a complete set of compatible operators J^, Jz, signyl and H, respectively. 
The dynamical operator A = a ■ {i + egr/r) — 1/2 is related to the Casimir of 0SP{1, 1), and 
the eigenvalues a = ±1 of signA describe the two helicity states of the electron related by the 
super conformal transformation. Hence the state of the system is given by J] 



H\j, m, a, E) = E\j, m, a, E) 



1 



1 



J \j,m,a,E) = j{j + l)\j,m,a,E), j = eg - -,eg + -,■ ■ ■ 
Jz\j,m,a,E) = m\j,m,a,E), m = + 1, - • • ,j - 1, j 



A\j,m,a, E) = adj\j,m, a, E), a = ±1, dj 



2 2 

e 9 
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The wave function in spherical coordinate and the Pauli two-component representation is 
'^E{r,0,4>) = {r,9,(l),a\j,m,a,E) = ^E{r)ilj>m,aiO,(p) Ul The angular part ??j,m,a(6', 0) can be 
expressed explicitly in terms of the monopole harmonics [3|] . 

Once the angular part of the wave function has been fixed, the Hamiltonian ([1]) reduces to 



1 1 d2 1 

H = — H ^ 

2m r dr'' 2mr^ 



i-A)i-A + l). 



(3) 



The radial eigenfunction ^Ei^) is the solution to the eigenvalue equation H<^e{i') = E(^e{t). 

We are interested only in bound states {E < 0) since the goal is to describe a system which can 
be used for quantum computing. The bound state radial eigenfunction reads 



<l>E{r) = Nr-^/^K2S,,^-i{f3r), 



P = {-2mEfl\ 5j,a = l-\a + ^d, 



(4) 



where Ki, is the modified Bessel function of the second kind and is a normalization constant. 
It is easy to see that only when u < 1, is normalizable over the region containing the origin. 
It turns out that there exists only one such bound state, ^E{f) = 2/3-fCi/2(/5^)/\/^) which arises 
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when j = \eg\ — 1/2. In particular, the orbit angular momentum is diagonal in the basis 
\j = \eg\ — 1/2, m, a, E) and there exists a ■ r/r\j = \eg\ — 1/2, m) = ziz\j = \eg\ — 1/2, m) J]. 

However, in this case ^E{f) is singular at the origin, and we need to perform regularization 
and renormalization operations on the Hamiltonian ([1]) to make it regular. It should be noted 
that the regularization of the model ^ is considerably complicated. First, the solvability of the 
theory depends on the larger dynamical symmetry 0SP{1, 1), which should be preserved as much 
as possible by the regularization scheme. Second, the theory has a U{1) gauge symmetry encoded 
in the angular wave function. The regularization should keep the angular part 'r]j.m,aid,(p) intact 
so that the Hamiltonian can reduce to the form ([3|). Otherwise, the exact solvability of theory will 
be ruined. The regularization scheme we use is described as follows. First, we observe that the 
reduced Hamiltonian ([3]) at j = |eg'| — 1/2 becomes 

1 1 (f 1 — e^g^ — ega ■ 
H = —prr + 



2m r dr^ 2m 

= -— V^-^fl + — U (5) 
2m 2m \ \eg\ J r'^ 

The Hamiltonian ([5]) implies that at j = \eg\ — 1/2 the radial sector of the spinning particle is 

equivalently described by a spinless particle in a spherically symmetric potential V{r) = — A/r^, 

A = (l + l/|e5|)eV/(2m). 

We choose a real-space cut-off regularization by introducing a length scale L as the regulator 

and re-defining the effective potential as Vji{r) = —X/r'^9{r — L), 9 denoting the Heaviside function. 

The regularized energy eigenvalue equation H(^e{t) = E(^E{f) reads 

" 2d \eg\i\eg\ + 1) 2mA 



'-^'v'-^' ' + —le{r -L)+ 2mE 



dr'^ r dr r'^ ' ^,2 \ ' • e{ ) ( ) 

The normalizable bound state solution expressed in the modified Bessel functions is 

^E{r) = ^r-i/2/|,^|+i/2(/5r), r < L, 

= Br-^/^K„^{l3r), r>L (7) 

The continuity of ^e{''') at r = L yields A = BKi^2{P^)/^\eg\+i/2{P^)j the normalization 

condition / drr'^\(f)E{^)\'^ = 1 fixes B = 2/3/y^ as L ^ 0. Finally, the continuity of dipEii") / dr 
Jo 

at r = L leads to 

Further, the expansion of Iu{x) near x = gives the lowest order reduction of ([8]) at L ^ 0, 
/3L = — {\eg\ -\- 1). Hence we get the regularized bound state energy E = -{\eg\ + lf/{2mL'^). 
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Obviously, the regulator dependent E is divergent as L ^ 0, and the spectrum is unbounded 
from below in this limit. There are two ways to cure this pathology. The first one is to adopt 
the viewpoint of the Wilsonian effective field theory We directly take the regulator L as the 
cut-off length scale A and consider the regularized Hamiltonian H = — V^/ (2m) + as an effective 
Hamiltonian above the length scale A. The bound state energy at r = A is 

The second one is the traditional approach of calculating the one-particle irreducible (IPI) effective 
action and performing a renormalization procedure as advocated in Ref. [7!]. At the renormalization 
scale r = A, we make the subtraction by splitting E = Eb + -Ediv In order to enforce the 
physical requirement that the wave function should vanish at the origin, which is needed for 
the self-adjointness of the Hamiltonian, we introduce a counterterm to the potential. This 
counterterm cancels the short-distance divergence 

Ediv = (l/A^ - (le^l + if /{2m) in the 

regularized energy E. Furthermore, as in field theory, the counterterm should be absorbed into 
the redefinition of the coupling constant A. One particular challenge in the present context is that 
the condition 2eg G Z must be preserved for quantum mechanical consistency. A detailed analysis 
of this procedure will be presented elsewhere j^. 

It is clear that both of the above approaches break the S0{2, 1) conformal symmetry due to the 
unavoidable presence of a length scale A. This is a direct manifestation of the conformal anomaly 
in this system, which has been shown in the modification of the SO{2, 1) commutator algebra 
through a deformation of the Hamiltonian by the anomaly operator: 

H^H + A, A = -i[H,D]A (10) 



The resulting anomalous conformal algebra 



BE 



9|] is composed of [H, D] = iH + [H, D]a, [H, K] 



2iD + 2t[H,D]ji and = iK + t'^[H,D]A. The Heisenberg equation further reveals the 

conformal anomaly as the non-conservation of the conformal charges, dD/dt = dK /dt = 2tA. 

An algebraic calculation of the first anomalous commutator shows that the anomaly operator 
A is directly related to the scaling behavior of the potential at the quantum level 7|, 

A = -i[H, D]a = i[H, D] + H=[l + h-V^ V{r) (11) 

We use the regularized wave function ([7]) and the regularized potential Vr to explicitly evaluate 
expectation value of the anomaly operator (jll|) . 

A = {A) = (y,s{T)) + \{r-VVR{r)) 
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drr^ ( 1 + -r— ) VR(r)|$E(r 



A straightforward calculation gives 



A 



m 



1 



1 + ^ ) = -2EBe^g^ ( 1 + ^ 



1 



(12) 



(13) 



On the other hand, an alternative and elegant interpretation on the origin of the anomaly in 
the Hamiltonian formalism has been presented in , where it was demonstrated that the anomaly 
is due to the fact that the symmetry generator does not leave the domain of definition of the 
lamiltonian invariant. By a careful observation on the Heisenberg equation, it had been shown 
} that the anomaly arises as A = i(^(t)| (h^ -H) G\^{t)), G denoting a certain symmetry 
generator operator which is D for the scale symmetry. This means that the anomaly operator can 
formally written as 

A={h^ -H)G = -iA\^{t)){^{t)\ (14) 

According to the rigorous definition of a self-adjoint operator [lo| . Eq. (I14j) implies that the Hamil- 
tonian has always acquired a non-self-adjoint piece once its domain of definition cannot be preserved 
by the symmetry transformation. 

The non-self-adjointness induced by the anomaly greatly modifies the quantum dynamics of the 
system. In the Heisenberg picture, the generator G satisfies a generalized Heisenberg equation Q], 



dG dG 



— = — + i[H,G]+i[H^-H)G= — +i[H,G]+tA 



dt dt 

and it implies the following time-evolution of G, 



dG 



G{t) = exp 



I j\s(H\s)- H{s)) exp iJ^dsH{s) 



G{0) exp 



dsH{s) 



(15) 



(16) 



In the Schrodinger picture, we have the time-evolution in terms of the modified Hamiltonian shown 
in Eq. ([lO]) 



\^{t)) = exp 



dsH{s) 



dsA\^{s)){'^{s)\ 



1^(0)) 



(17) 



The formal integration solution (|17|) for |^'(t)) shows that in the presence of the anomaly the 
quantum system undergoes a non-unitary evolution resultant from the anomaly. This is consistent 
with the fact that anomalous effects in a quantum theory contribute only to the imaginary part of 
the quantum effective action |ll| . 

Turning to the model at hand, we take G to be the generator D of the scale symmetry. The 
conformal anomaly arises only for the normalizable bound state ^'£;(r, 6, cj), a) in the s-wave sector. 
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and originates from its radial part $£;(r). Therefore, Eq. (jl7p tells that the time-evolution for this 
specific stationary state should be 

\^E{t)) = e-*(^-*-^)*|<I>i=;(0)) (18) 



where the energy and anomaly are provided by Eqs. and (jl3p . respectively. Note that the 
completeness condition = 1 is used in deriving Eq. ([18]). Although all the energy 

E 

eigenstates, including the scattering states, should be taken into account in the the completeness 
condition, the anomaly only pertains to the bound state and vanishes for all other eigenstates. 
Thus we effectively take \^e){^e\ = i- and get Eq. (fTHj) . 

We now consider the electron-monopole system as a physical model for quantum computing. The 
quantum state we manipulate is just the normalizable bound state ^£;(r, ^, a) at j = \eg\ — 1/2, 
its two-level spin degrees of freedom playing the role of a qubit: 




The spatial amplitudes fi{r,6,(j)) {i = 1,2) can be obtained with some algebraic operations J]. 

As we will show below, it is the time-evolution of fi{r,6,(j)) related to ^e{t) that brings about 
the decoherence between two spin states during a quantum computing due to the presence of 
anomaly. Roughly speaking, the two spin states constitutes a qubit, and one must control their 
dynamical evolution to carry out information processing. We therefore switch on a time-dependent 
Hamiltonian to make the spin flips that can ultimately be used in a quantum algorithm. However, 
the spatial sector fi{r,6,(f)) of the wave function will evolves in time controlled by the quantum 
effective Hamiltonian of the system itself along with the spin flipping dominated by the external 
Hamiltonian. According to Eq. (jl8p . the anomaly will cause /j(r, 0, (/>, t) to have a damping factor 
which in turn will lead to decoherence. In the following we show the details of how this phenomenon 
happens. 

Let us first analyze the quantum effective Hamiltonian provided by the system itself. Obviously, 
the time-evolution ([TH]) of ^e{t) of the bound state wave function gives the spatial part, -ffgpa = 
E — iA. As for the spin sector, we use the fact that at j = \eg\ — 1/2 the orbit- and spin- angular 
momenta decouple, and the spin part of the wave function is the eigenstate of the operator a ■ r/r. 
Specifically, the form of the radial Hamiltonian ([5]) shapes only when the eigenvalue equation of 
the operator a ■ r/r has been applied. So we can simply choose -ffgpin = f ■ ^/^- A combination 
of the spatial and spin sectors determines that the effective Hamiltonian with resect to the bound 
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state ()19p should take the following form: 



sys 



spm 



{E - iA) 



a ■ r 



{E — iA) {ax sin i9 cos ip + ay sin i9 sin ip + az cos i?) 



(20) 



where ('i?,'/?) represents the spin orientation in three-dimensional space. 

Eq. (j20p is the effective Hamiltonian realized on the bound sate of the system. We now switch 
on a time-dependent external Hamiltonian to make the spin flip. A typical choice is the inter- 
action of the spin with an oscillating external magnetic field in two-dimensional x — y plane, 
Bext = Bq (coswte^. + s'mujtex), and the Hamiltonian -ffcxt = e/2'ma ■ Bcxt- The spin dynamics is 
dominated by the Schrodinger equation 



' dt 



'eB, 



eBo 
2m 



cos ijjt + {E — iA) sin -d cos ip 



a^. 



+ 



— ^sinwt + - i^) sin??sin99 a^^ (E - iA)co's.'Qaz\^Eit) (21) 
2m J J 



We neglect the E — iA term in the ax and ay components since usually the microscopic values of 
the energy E and the anomaly are much smaller than the macroscopic magnetic field, |^| <C 
|e|-Bo/2m. In this approximation the time-evolution of the spin state reads 



^E^S) = exp < — Acq's, -da zt 



UJ 



E cos '& — — ] az + 



eBo 
2m 



t Ub(0) 



(22) 



To show explicitly the occurrence of the decoherence implied from ^ E{t), we take Ecos'd = (jj/2 
as in nuclear magnetic resonance and use again |^| <C \e\BQ/2m. Assume that the initial state is 
spin-up, ^'s(O) = /i(r, e, 0) (1, 0)^, Eq. ^ yields 



■^E{t) = Ci{t) 



+ C2{t) 




Clit) 



cos 



^-A'^cosH] ' t 
4m^ / 



sm 



A cos ?? ■ 



4m^ 



1/2 



sm 



C2(t) 



2m 



^-A'cos'^^ 



1/2 



4m ^ 



1/2 



(23) 



Clearly, the non- vanishing A leads to |ci(t)p + |c2(t)P / 1, and hence the decoherence between 
the two helicity states occurs and the qubit is destroyed. 
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To summarize, we have used an electron-monopole system to reveal a phenomenon not previ- 
ously discussed in the quantum computing literature: a quantum mechanical anomaly can result 
in decoherence. Note that anomaly is a quantum dynamical phenomenon rooted within the sys- 
tem itself. It reflects how quantum efl'ects can render a classically feasible symmetry unrealizable. 
One typical example is the case where the configuration space has non-trivial topology so that 
the Hilbert space constructed via the quantization procedure from the classical phase space can- 
not sustain all the classical symmetries. In the case we have just considered, the source of the 
anomaly is the singular behaviors of the interaction potential near the magnetic monopole. The 
classical conformal symmetry does not preserve the Hilbert space as the domain of definition of 
the Hamiltonian due to the singular behavior of the wave function in the s-wave sector. 

Until now the search for a physically realizable quantum computer has been concerned only 
with decoherence that arises due to interactions with the external environment. It is important 
to emphasize that decoherence can also in principle be induced by quantum anomalies. Since this 
dissipation originates from the dynamics of the quantum system itself, it seems that it has the 
potential of being more destructive than the standard mechanisms for decoherence. 

This work is supported by the Natural Sciences and Engineering Research Council of Canada. 
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